Methods of optimizing X-ray optical prescriptions for 

wide-field applications 

Ronald F. Eisner", Stephen L. O'Dell", Brian D. Ramsey", and Martin C. WeisskopP 
■NASA Marshall Space Fhght Center, Space Science Office, VP62, Huntsville, AL 35812 



ABSTRACT 

' We are working on the development of a method for optimizing wide-field X-ray telescope mirror prescriptions, 
, including polynomial coefficients, mirror shell relative displacements, and (assuming 4 focal plane detectors) 
' detector placement along the optical axis and detector tilt. With our methods, we hope to reduce number 
^ , of Monte-Carlo ray traces required to search the multi-dimensional design parameter space, and to lessen the 
complexity of finding the optimum design parameters in that space. Regarding higher order polynomial terms 
I as small perturbations of an underlying Wolter I optic design, we begin by using the results of Monte-Carlo ray 
' traces to devise trial analytic functions, for an individual Wolter I mirror shell, that can be used to represent 
the spatial resolution on an arbitrary focal surface. We then introduce a notation and tools for Monte-Carlo ray 
tracing of a polynomial mirror shell prescription which permits the polynomial coefficients to remain symbolic. 
In principle, given a set of parameters defining the underlying Wolter I optics, a single set of Monte-Carlo ray 
^ ■ traces are then sufficient to determine the polymonial coefficients through the solution of a large set of linear 
equations in the symbolic coefficients. We describe the present status of this development effort. 

Q ■ Keywords: X-ray astronomy. X-ray optics, ray trace, wide field-of-view optimization 

S ; 1. INTRODUCTION 

In 1992, Burrows, Burg, and GiacconP' showed how by adding higher order polynomial terms to Wolter I 
prescriptions, and hence giving up some on-axis spatial resolution, one can obtain prescriptions for reflecting 
surfaces that provide improved average spatial resolution over a wide field-of-view (say ~ 30 arcmin). Such 
so-called polynomial optics would be particularly useful for moderately deep to deep surveys, to be carried out 
! by observatories such as for the proposed Wide-Field X-ray Telescope (WFXT) mission,'^! and for solar X-ray 
l/^ ■ observations. Procedures for optimizing the design of wide-field X-ray telescopes utilize Monte-Carlo methods 
for determining the design parameters, including specification of the polynomial coefficients.'^'^'^ Monte-Carlo 
' ray traces are performed over a range of design parameters, and the final design determined according to some 
, optimization criterion and methods. Since the number of mirror shells per module is typically large (~ 50 — 100), 
these procedures are presently complicated and computer intensive. 

The present paper is a report on the current status of an on-going studj^ of the properties of Wolter I and 
, polynomial optical prescriptions with the ultimate goal of simplifying the procedures for optimizing their designs. 
^ ' Since a polynomial prescription can typically be viewed as a small pertubation to an underlying Wolter I design, 
5^ , we begin with Monte-Carlo studies of the properties of Wolter I mirror shells relevant to wide-field designs, 
attempting to deduce analytic formulae for representing the geometric area and spatial resolution as functions 
of source position on the sky relative to the pointing axis, focal length, mirror shell segment length and shell 
intersection radius. We then outline a method, valid when the polynomial coefficients are sufficiently small, for 
ray tracing polynomial optics keeping the polynomial coefficients in symbolic form. 

A merit function providing a measure of spatial resolution averaged over the field-of-view (FOV) is defined 
in f|21 while the parameter space over which we have carried out Monte-Carlo ray traces is described in ^ 
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We note in ^ll that the spatial resolution, when averaged over the FOV as in the merit function, as a function 
of source position relative to the optical axis, is a simple sum of terms up to second order in (1) the mirror 
shell displacement relative to the nominal on-axis focus, and (2) the tilt angle for the CCD detector array. In 
[JS] we arrive at the important conclusion that the spatial resolution on an arbitrary focal surface for a set of 
nested mirror shells may be written as the sum of two terms. The first is a sum over the spatial resolution 
of the individual shells on that surface, weighted by their effective area. The second is a sum over a kind of 
weighted variance of the mean ray positions for the individual telescopes on that surface. In ^ we introduce a 
compact notation for representing ray trace variables such as position or direction vectors, including polynomial 
coefficients in symbolic form. In ^ we discuss the outer product of two vectors, a concept from linear algebra 
necessary for the development of the methods introduced in this paper. In ^ we specify the basic operations 
of a polynomial optic algebra, which are addition, subtraction, multiplication, division, and the taking of square 
roots. Given a direction vector, ki, and initial position xi = {xi,yi, zi), ^j9]shows how to propagate a ray from 
sci to axial position Z2 and determine the other coordinates, X2 and y2, thus determining the final position a?2. 
We define our coordinate system and the mirror surface prescriptions for polynomial X-ray optics in ^101 In 
we list the tasks required to trace rays through X-ray optics. In the future, we plan to show how the tools 
presented in this paper are used to accomplish these tasks while keeping the polymonial coefficients in symbolic 
form, and to provide concrete examples. Some closing remarks are provided in 



2. MERIT FUNCTION 

For X-ray survey applications, such as the proposed Wide-Field X-ray Telescope (WFXT) mission,'^! one desires 
a large effective collecting area over a broad energy range combined with good spatial resolution over a wide 
FOV. The geometric area available is essentially pre-determined by the diameter of the launch vehicle faring, 
the number of desired telescope modules (which are constrained by the desired FOV and, in the absence of 
extendable optical benches, the focal lengths permitted by the launch vehicle faring), and the number of mirror 
shells per module allowed by mass and manufacturing constraints. In our work, we have therefore concentrated 
on optimizing the spatial resolution average over the FOV, by minimizing the merit function: 

where 9 is the polar off-axis angle for the incident X-rays, is the azimuthal angle for the incident X-rays, and 
w{9, (j)) is a weighting factor. By symmetry, the average in Eq. ([T]) may be restricted to ^ S [0,7r/4] for a typical 
detector setup consisting of four tilted CCDs, each occupying a single quadrant. This statement neglects any 
repositioning of the detectors to place the on-axis aim point on one of them. The quantity (7^(0, 0) is the variance 
in the position of rays reaching the focal surface. This focal surface may be curved or tilted with respect to the 
flat plane perpendicular to the optical axis and passing through the nominal on-axis best focus. The variance, 
a^{9, (j)), is given by 



^2(61, 0) = [ ( < a;2 > - < X >2 ) + ( < y2 > _ < y >^ ) + { < > - < z >^ ) ] (2) 

where x, y and z are the positions of the rays on the chosen focal surface, and < q > denotes an average of the 
quantity q. All rays incident on the detector are included in the averages in Eq. independent of the mirror 
shell from which they exited. 

We have found that the coefficients of the polynomial terms modifying Wolter I optics for wide-field appli- 
cations may be regarded as small for our purposes. Therefore, we treat them as small perturbations to the 
underlying Wolter I design. It is for this reason that we have sought analytical fitting functions for the contri- 
butions to (j'^{9,(p) for Wolter I optics. This is also the justification for the procedure we describe later in this 
paper for ray tracing polynomial X-ray optics keeping the polynomial coefficients as symbolic and unevaluated 
until optimized. 




Figure 1. Mirror segment length Is vs. intersection radius ro,s, with points showing locations in the (rj ,,,^3) plane 
of Monte-Carlo ray traces with 50,000 incident rays for focal lengths of 5.5 m. The largest dot shows the location of 
additional ray traces with 100,000 incident rays at a focal length of 5.5 m. The largest and mid-size dots show locations 
of additional ray traces with 50,000 incident rays for focal lengths of 4.5, 5.0 and 6.0 m. Note there are 2 segments 
per shell. The solid line represents our reconstruction of the wide-field telescope design described in Ref. [7| using their 
design constraints. The vertical dashed lines show constant values for nominal graze angles of 25, 50 and 90 arcmin at 
the intersection plane for chosen values of ro. The curved dashed lines show constant values of 5 and 15 arcmin for Ocoma 
[see Eq. (fT^ ]. in the (4,^0,3) plane. 

3. MONTE-CARLO RAY TRACES 

In order to explore the dependences of <J^{9, <p) on mirror shell parameters, we carried out an extensive series of 
Monte-Carlo ray traces of single shell Wolter I optics for nominal focal lengths, /, of 5.5 m, mirror shell segment 
lengths (2 segments per shell), £s, of 10, 15, 20 and 40 cm, and shell intersection radii, rp^s, of 15, 30, 45 and 60 
cm. Here the subscript s denotes a shell number. Figure ([1]) plots the locations of these ray traces in the £s vs. 
ro,s plane. In most cases, the number of rays incident on the shell aperture was 50,000; for the point in Figure 
^ marked with the biggest dot the number of incident rays was 100,000. The mid-size dots show locations of 
ray traces for focal lengths of 4.5, 5.0 and 6.0 m. We used the results for from these ray traces to devise trial 
analytic functions for representing cr^(0, 4>) a-s a function of the angles 9, and of /, is and tq^s. 

The solid curve in Figure ([T]) shows our reconstruction of the relationship for is vs. tq^s for the 3 telescope 
module, 82 mirror shell per module wide-field design described and discussed in Ref. [T] We carried out this 
reconstruction using the design constraints provided in their Table (2). While certain of our assumptions may 
vary from theirs, in general we expect our reconstruction to be close to their actual design. 

4. SINGLE MIRROR SHELL 

For Monte-Carlo ray traces of a single mirror shell s, we define the geometric area, Ageom.s, as 

-Ageom,s{9) = -^inc^s / 1T^inc,s{9\ (3) 



where Ainc,s is the entrance aperture for shell s, nines the number of rays incident on that aperture, and Us the 
number of doubly reflected rays exiting the mirror shell. 

On a detector tilted by an angle 0tiit with one corner at the {x, y) origin, but displaced along the optical 
axis by an amount 5zs from the flat plane perpendicular to the optical axis at the nominal on-axis focus, the 
variance, or square of the RMS dispersion may be written in the form 



a'^{0,(j),Szs,0tiit) = cis +2bsSzs + Cs Sz"^ + 2 tanBuit + 2 Cg Szg ta.ii9uit + fs t&v? Quit- (4) 
Evaluating the merit function [Eq. ([T])] for this shell leads to 

M{6zs,9tiit) as^M + 2bs.M6zs^M + c^.m Sz^;^ 

+ 2 4,M tan 6*41/4 + 2 Cs^m Szg lanOtat + fsM ta.T? Otm, (5) 

where the subscript M denotes an average over the FOV like that in Eq. ([T} . In order to carry out the integrals 
over the FOV, it is advantageous to have analytic forms for the coefficients as, bg, c^, dg, and fs- Minimizing 
Eq. Q with respect to dzg and tSLnOtut, we find 

n I i's,M &s,M — Cs^M ds ]\I \ 

tan^tiit = 2 — (6) 

Y Cs,M Js,M - e^ j^j J 
r. I ds^M es,M — bs^AI fs.M \ 

SZs = 2 — (7) 

y Cs,M Js.M - e^ j^,j j 

Expressions (|4]) — ([7]) are general and applicable to any surface prescription for grazing incidence X-ray optics. 
Below we use the notation: 

5:<,„,(^|.)>. <,(y>, 

and similarly for other combinations of terms. Angle brackets around a quantity, < q >, denote an average over 
that quantity on an the focal surface. Angle brackets with the subscript 0, < q >o, denote an average in the flat 
plane perpendicular to the optical axis at the nominal on-axis focal position. We find that the coefficients as, 
bs, Cs, ds, Cs and fs can be expressed in terms of averages in that fiat plane. Assuming a detector in the first 
quadrant (both x and y positive) , here we provide some examples of coefficient definitions along with the trial 
fitting functions that we find useful for representing the results of our Montc-Carlo ray traces: 

^2(0,0) = as = i<i^:yf>o,s - <ix,y)>l, ), (9) 

The behavior of Us as a function of 6 is complicated by effects due to coma. At present, we are working with the 
three trial fitting functions, ao{9), ai{9) and 02 (6*), for Wolter I optics. We define a useful function, g, and then 
ao{9), ai{9) and 02(6*): 



9{d, C, C) = 1 + C tane + C ta.n^9 



(10) 



afitiO) = ttcomaiO) + ttmiO), (to = 0, 1, 2) 



aco,na{0) = (tan4ao / 2 tan^ 9 aQ{0) = ( 2 ^a,o ^ / tan4ao f tan^ 6* g{e, (a, Ca) 

ai{9) = acoma{9) + ao{S) 02(6') = (^V ^comaiO) + Vo^jM) > 



(11) 



We note that ai defined here is equivalent to Eq. (2) in Ref. i9i Additional examples of coefficient definitions 
and trial fitting functions applicable to Wolter I optics are: 



>o,s - < ix,y) >o,s< 



>o„ 



(12) 



E 



< 



>0„s - < 



^{x,y) \ 2 



CfuiO) 



tan4ao g^O, (c Cc) f 



(13) 



Definitions and trial fitting functions for ds, , and fs are lengthy, so for reasons of readability we provide them 
in Appendix 1X1 

5. NESTED MIRROR SHELLS 

Consider a set of S nested telescopes. We assume uniform illumination of the entrance aperture for the full array 
of telescopes. The number of rays through the s-th telescope is n^, and the total number of rays through the 
array of nested telescopes is 



N 



E^ 



(14) 



We designate the (x, y, z) position coordinates on an arbitrary focal surface for the full array of the k-th ray 
through the s-th telescope by ( x^^^, z^^k ). We find that the total variance, or square of the RMS disperion, 

for the full set of nested shells may be written in the form 



^2 , 2 



(15) 



where 



(16) 



with a1{5zs, Ouu) given by Eq. (g]). The second term on the right-hand-side of Eq. p5)) is given by 



N 



N -1 



E 



S 

'n 

s—1 



/ s 



(17) 



Eqs. ([15]) — (|17l) show that the variance, cr^, for the full set of nested shells has two contributions. The 
first, Eq. (jl6p . is a sum over the variances for the individual telescopes, on the chosen focal surface, weighted 
essentially by their relative effective geometric areas [(n^/iV) ~ ^geom.sl^*)/ ^9eom,s(^)]- The second, Eq. 
(jl7p . is a sum over a kind of weighted variance of the means, < (a;, y, z)s >, for the individual telescopes on that 
focal surface. This second contribution can be viewed as arising from the differences in the best focal surfaces 
for the individual mirror shells from that for the full set of nested shells (see Ref. [7]). Expressions — (|17l) 
are general and applicable to any surface prescription for grazing incidence X-ray optics. 

For best performance, Eqs. (|15p — (jl7p mean that the minimization of cr^j, and thus the optimization of the 
parameters Otm and {Szs,s = 1, 2, 3, ...S*), must be done simultaneously, rather than following Eqs. ^ and ([7]) 
for the individual shells. In principle this can be done using matrix methods, although the number of linear 
equations involved is large for current wide field designs which approach 100 nested mirror shells (see Ref. [7]). 

We have also shown a need for expressions for < {x,y,z)s >, and terms, < {x,y,z)s >^ and cross terms 
< {x,y,z) X {x,y,z) >, that can then be derived to the appropriate order, for the individual shells. We find 
to the appropriate order 



(18) 



< z >s = ( a^,s + a'y ,, )taiieuit + h'^ .. + ( by ,, ) 5zs i&nduu + ( d'^ ,, + d'y ,, ) tan^ 



We provide the definitions of a[j.^y)^s, ^[x,y),s^ '^'{x,y).s^ '^'{x.y).s f[x,y),s AppendixlX] For use below, we define 



an angle, Ocoma, at which acoma{0) and ao(0) [see Eq. pT|) ]. with Ca = and ~ 0, are equal: 

tan 6'coma = (^) (l) t^^^^ao. 



(19) 



We have not yet finished devising analytic expressions for e'^^ y^ ^ and for /^'^ ^, but here provide trial fitting 



functions for for a', % , n and d', s . 

(a:,j/),s " {x,y),s 



"■'(x,y)Jit = / ( 1 + 5f(x,y),a ) tau ( 



tan^ e 



5(^5 C{x.,y) .ai ^{x.y) .a) 1 



(20) 



{x,y)Jit,OV'^' 



land (cos0, sin0) g{9, C(x,y),b, ^{x,y)M) 



^(x,y),coma{^) — P(x,y),coma sin [ TT C{x,y),coma {P I ^ coma) 
^{x,y),damp{(^) ~ GXp [ ^{x,y),coma ((^ /Scoma) ] 



(21) 



^[x,y)Jiti^T 4') — fc(2;,a),(iamp(^') ^(a:,a),coma(^) + (1 ~ ^(2;,y),damp('^)) [Q{x,y),coma + ^[x ,y) Jit, oi^ ^ 'i') 



d', 



(x,y),ntiS,^) = - / (1 + '5/(x,,),rf)(l + 2 tanMao) tan^ 



(22) 



X (cos(/), sincf)) {coscj) + sincf)) g{e, C.{x,y),d^ k^,y),d)^ 
In the future, we plan to provide a fuller account of our methods and results for ray tracing Wolter I optics. 



6. NOTATION FOR POLYNOMIAL COEFFICIENTS 

In the past, we have written sums over rays 



A - ^Afe, (23) 



fe=i 



where A^ is some quantity such as position along an axis or a component of a direction vector for ray k, in the 
form of a second order expansion in the polynomial coefficients 



A — Aqooo + ui Aiooo + U2 Aoioo + U3 Aooio + U4 Aoooi 
+ Ui A2000 + U2 A0200 + A0020 + u1 A0002 

(24) 

+ ui U2 Aiioo + Ml M3 Alolo + ui Ui Aiooi 

+ U2 U3 Aoiio + U2 Ui Aoioi + M3 Ui Aooii, 

where Ui, U2, W3, and U4 are the polynomial coefficients for the mirror shell. The form Eq. (I24p assumes the 
coefficients are small enough so that a second order expansion is valid. The polynomial deviations from Wolter 
I optics required in the applications we have studied so far satisfy this criterion. 

We define a polynomial coefficient vector 



u = {ui,U2,U3,Ui). (25) 



We also define the scalars 



Aqo = Aoooo Aqi = Aiooo A02 = Aqioo 
Ao3 = Aooio Ao4 = Aoooi An = A2000 

A22 = A0200 A33 = A0030 A44 = A0002 (26) 



2 Aiioo A13 — A31 — ^Aioio Ai4 — A41 — 



A23 — A32 = -Aoiio A24 — A42 — ^Aoioi A34 — A43 = -Aooii, 



and the vectors 



Ao = (Aoi, A02, Ao3, A04) 
Ai = (All, A12, A13, A14) 



A2 = (A21 , A22 , A23 , A24) 

As = (A31 , A32 , A33 , A34) 
A4 = (A41, A42, A43, A44). 
Finally we define a matrix with the row vectors Ai, A2, A3, and A4: 



A 



(27) 



( X^] 




1 All 


A12 


Al3 Ai4 


\ 






A21 


A22 


A23 A24 


h 




A31 


A32 


A33 A34 




\ A4 J 




V A41 


A42 


A43 A44 


/ 



A = Aoo + w • Ao + u ■ X - u. 



(28) 



(29) 



In this notation, and since the matrix A is symmetric, derivatives of A with respect to tlie polynomial coefficients 
are 



dX 



Aoj + 2 ^ Xji Ui = Xoj + 2 u - Xj . 

i=l 



For a set of nested shells, denoted by prefixes s or r {e.g., Xg), we note that 

dX. 



d A^^ /dusj 



dusj 



7. OUTER PRODUCT OF TWO VECTORS 



(30) 



(31) 



Consider the two vectors 



a = (ai, 02,03, 04) 

b = {bi,b2,b3,b4). 
In linear algebra, the outer product of b with a is given by 



6 (8) a 



( b\ ai b\ a-i b\ b\ 04 \ 

62 a\ 62 0,1 62 03 62 04 

63 oi &3 ai 63 03 63 04 
y 64 a\ &4 ai 64 a-i 64 04 ) 



The outer product is an essential tool in polynomial ray tracing algebra. 



(32) 



(33) 



8. BASIC OPERATIONS 

Consider the two polynomial objects 



a = aoo + u- do + u-a-u 

b = boo + u ■ bo + u -b ■ u (34) 

U = {ui,U2,U3,U4). 

Treating il as small and expanding to second order in il, we now specify how to carry out basic operations 
on a and b. Using these operations, it is possible to construct a Monte-Carlo ray trace code for 
polynomieil X-ray optics with sufficiently small but unknown coefficients u. In principle, values 
for the coefficients can then be derived from a final ray bundle for any assumed merit function. 

In the case of the addition operation, we note 

a ± b = {aoo ± ^'oo) + m • (oq ± 60) + u ■ {a±b) ■ u. (35) 

The multiplication operation is more complicated. We want to keep terms only to 2nd order in the polynomial 
coefficients u. To this order we find 

axb = aoo &00 + u- (aoo ^0 + ^'oo oo) + u- (aoo 6 + 600 i + &o ® ao) • u. (36) 
In particular, we note that 

{il ■ do) X {u- bo) = u- {bo <E) do) ■ u. (37) 
We also need the square and the cube of a polynomial object to 2nd order in u. We find 

«^ = oioo + ^ ■ (2 aoo oo) + w • (2 aoo a + ah® ah) ■ u 

(38) 

a^ = ^00 + ^ • (3 "oo dh) + u - (3 ago a + 3 aoo dh (8) dh) ■ u. 
For the division operation, to second order in u, we find 




(39) 



In the case of the square root operation, to second order in u, we find 



We also find 



a I ^^h = 



300 



flo 
floo 



boo 



«oo 



000 



- 

V aoo 



&00 



boo 



ao 
aoo 



b^ 

boo 



(42) 



b^ 

boo 



Eq. (j42p is especially useful for normalizing unit vectors such as ray direction vectors or normal vectors to 
surfaces. 

9. PROPAGATION OF RAYS 

The new x and y coordinates of rays propagated from axial position zi to Z2 according to 

X2 = Xi + k^i t 



y2 



yi + kyi t 



(43) 



Z2 



zi + kzi t. 



Since zi and kzi are known, the parametric variable t can be expressed in terms of those quantities and the value 
for Z2. In polynomial notation, we have 



t — too + u ■ t + u ■ t ■ u 



Z2 - Zi 

kzi 



[ (^2,00 - 2i,oo) + u-{z2- zi) + u ■ [z2 - zi) ■ u] I {kzx^oo + u ■ kzi + u ■ kzi ■ uj . 



(44) 



The rule for division, Eq. p9l) . shows how to evaluate the polynomial components of t namely (too, t, and t) 
from Eq. ()44|) given the known polynomial components of zi, Z2, and ki. Then, also given xi, yi, k^i, and 
kyi in the form of polynomial objects, we can apply the basic operations of addition and multiplication defined 
in Sj8] to Eq. (|43|) to find X2 and 2/2 in the form of polynomial objects also. This means that knowing the 
polynomial components of the initial position and direction vector, we can compute the polynomial 
components of the final position vector without knowing numerical values for the polynomial 
coefficients u. 



10. SURFACE PRESCRIPTIONS FOR POLYNOMIAL X-RAY OPTICS 

We consider mirror prescriptions for the primary mirrors of the form 



= rl^ l + 2As (z/ro^s) + {z/ra^sf + (^/ro,.)^ • Co!, 
For the primary (P) and secondary (S) mirror segments, we have 

P : ^=tanao,s, 5. = 0, Co,s,p = (l, (7^) , 0, 

5*: As=tan3ao,s, -B^ = /i(ao,s) tan^ 3ao,s, Co,s,s = (o, 0, 1, ^ 

Here 



ao,s = I ^ tan ^ 



(45) 



(46) 



(47) 



h{aa,s) = 1- l/[l + 2cos(2ao,s)]', 



and u — {ua^s,PTUijgp^Ua,s,s^'^b,s,s)- The notation and methods introduced here are, in principle, readily 
extended to additional terms in the mirror presciptions [e.g.^ proportional to (z/ro,s)^, [z/tq^s)^ , etc.]. 

11. RAY TRACING POLYNOMIAL X-RAY OPTICS 

In order to trace rays through X-ray optics, one needs to: (1) populate the entrances aperture with rays (both in 
position and direction); (2) calculate intersections with mirror segment surfaces; (3) calculate the unit normals 
to the surfaces at those intersections, including deviations due to non-ideal surfaces; (4) determine the direction 
of the reflected ray; and (5) take account of obstruction by the next innermost mirror shell. Using the tools 
outlined above, all these tasks can be accomplished, to suflicent accuracy, for polynomial X-ray optics as long as 
the polynomial coefficients are sufficiently small. In the future, we plan to provide more details on this method 
and its results. 



12. CONCLUDING REMARKS 

The ultimate goal of our work is to reduce the complexity of design procedures for nested grazing incidence X-ray 
telescopes, specifically those with Wolter I and polynomial designs. In this paper, we have: 

1. Described our use of Monte-Carlo ray traces to devise trial analytic formulae for the coefficients of terms in 
the expression [Eq. ^] for the spatial variance of rays from a point source on an arbitrary focal surface for 
a single Wolter I mirror shell. Our adopted merit function [Eq. ([1}] can then be minimized to provide the 
best displacement of the mirror shell along the optical axis and the best value for the detector tilt angle. 

2. Shown that for a set of nested mirror shells, the spatial variance on an arbitrary focal surface is a sum 
of two terms. The first [Eq. ([TE| ] is a sum over the variances of the individual shells evaluated on that 
focal surface, weighted by their relative effective areas. The second [Eq. (IT71) ] is a sum over a kind of 
variance for the mean positions of rays from the individual shells on that focal surface. The existence of 
this second term means that it is necessary to optimize parameters such as mirror shell displacement along 
the optical axis, detector tilt angle, and polynomial coefficients simultaneously for all mirror shells, rather 
than individually. 



3. In ^S] — ^111 introduced notation and mathematical tools for ray tracing polynomial optics leaving the 
polynomial coefficients in symbolic form. In principle, this simplifies the design procedure by reducing 
the required number of Monte-Carlo ray traces, and permitting determination of numerical values for the 
polynomial coefficients through the solution of a large number of linear equations derived from minimization 
of the merit function. 

Our future plans are to continue these studies, refining the trial analytic functions for Wolter I optics, 
implementing a polynomial optic ray trace code using the tools described in this paper, and hopefully providing 
a less complex means for the optimization of wide-field X-ray telescope designs. 
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APPENDIX A. ADDITIONAL COEFFICIENT EXPRESSIONS 

We assume a detector in the first quadrant (botli x and y positive) . First we provide definitions for tlie coefficients 
ds, Cs and fs (see g]): 



ds ^ < (^'2^) f%^) (^ + 2^) >M - < {^,y) >o,s< ^^^"'^'^ 



{x + y) >o,5 



(48) 



E 



< 



{x + y) >o,s - < 



>0,s< 



{x + y) >o„ 



(49) 
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<(:. + y)^(^)'>o,. - <(. + y)(^) 
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0,s 



2 E 



< ix,y){x + y) 



>o,s - <{x,y)>^,s<{x + y) 



"^{x.y) 



(51) 



kx ~\~ 



/^.^ =<{x + yf>o - <x + y>l 



We typically find f^^s ^ We define 



g{0, C, = 1 + C tan6' + ^ tan^ 



(52) 



ho{<j), rj) — 1+7] sin 20 



(53) 



HO,(l), C, 'yo, = ^o(0, ?/o) + C tan 61 /io(0, »yc) + ? tan^ 6* /io(0, ?/?)• 

Note that h{9,(p, ^, 0, 0, 0) = g{(^,(f>, 0- Now we provide trial fitting functions for the coefficients dg, 
and fg-. 



dfr. 



2 fid f i tan^e h{9,(t), Q, f^, ?7o,d, Vcd, Vi,d), 



(54) 



efit{9,(l)) = / tan^4ao i&nO h{e,(t), (e, ■fe, Jyo.e, nCe-, f?c,e), 



(55) 



fxyjit{0,(t>) = 2 ( / tan4ao tan^ f h{d,(t>, Cf,xy, ^f,xy, mj,xy, VCJ,xy, Vi,f,xy), (56) 



fzjit{0,4>) = /i/,^ tan^ 6* /i(6', (/), C/,^, Ce, Vc,f,z, ViJ,z)- 



(57) 



We now provide the definitions of a[^ y^^^, ^[x,y),s^ '^[x,y),s^ '^[x,y),s and /('^ ^ (see gS]): 



0'{x,y),s = < >0,s 



d[x,y),s =<(%f)(^ + 2/)>0,s, (58) 



